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£NJ . Abstract: In this paper we consider a viscoelastic wave equation with a time- 

varying delay term, the coefficient of which is not necessarily positive. By intro- 
ducing suitable energy and Lyapunov functionals, under suitable assumptions, we 
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establish a general energy decay result from which the exponential and polynomial 



< 

o 

■ types of decay are only special cases. 
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1 Introduction 



' In this work, we investigate the following viscoelastic wave equation with a linear damping 

^ ; 

, and a time- varying delay term in the internal feedback 

O 1 
OO 

o 



r ^ 

uu(x, t) — Au(x, t) + / g(t — s)Au(x, s) ds 
Jo 

+aoUt(x, t) + aiut(x, t — r(i)) = 0, (x, t) G x (0, oo), 
u{x, t) = 0, (x, t) G dn x [0, oo), ( L1 ) 

u(x, 0) = Uq(x), ut(x, 0) = u\(x), x&Q, 
{ u t (x,t) = f {x,t) (x,t) e!lx [-r(0),0), 

where $7 is a bounded domain of M n (n > 2) with a boundary <9Q of class C 2 , ao and a\ are 
real numbers with ao > 0, r(t) > represents the time-varying delay, and the initial datum 
uo,u\,fo are given functions belonging to suitable spaces. 

The viscoelastic wave equation without delay (i.e., a\ = 0), has been considered by many 
authors during the past decades. Cavalcanti et al. [7] studied 

u tt -Au+ / g(t — r)Au(r)dr + a(x)u t + |ii| 7 u = 0, (x,t) G O X (0, oo), 
J o 

for a : 17 — > M + , a function, which may be null on a part of the domain f2. Under the 
conditions that a{x) > ao > on uj C f2, with cj satisfying some geometry restrictions and 

-£i9(t) < g'(t) < -&g{t), t > 0, 
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the authors established an exponential rate of decay. Berrimi and Messaoudi [3j improved 
Cavalcanti's result by introducing a different functional which allowed to weak the conditions 
on both a and g. In [8J, Cavalcanti et al. considered 



under similar conditions on the relaxation function g and a(x) + b(x) > p > 0, for all x G 17. 
They improved the result of [7| by establishing exponential stability for g decaying expo- 
nentially and h linear and polynomial stability for g decaying polynomially and h nonlinear. 
Berrimi and Messaoudi [1] considered 



in a bounded domain. They showed, under weaker conditions than those in [8], that the 
solution is global and decay in a polynomial or exponential fashion when the initial data 
is small enough. Then Messaoudi [19j improved this result by establishing a general decay 
of energy which is similar to the relaxation function. For other related works, we refer the 
readers to [SI El EH ESJ \M EH H3 ED E3 E3] and the references therein. 

In recent years, the control of PDEs with time delay effects has become an active area of 
research, see for instance [T| 1161 [771 128] an d the references therein. The presence of delay may 
be a source of instability. For example, it was proved in [91 [TOl [2T1 [22"1 [30] that an arbitrarily 
small delay may destabilize a system which is uniformly asymptotically stable in the absence 
of delay unless additional conditions or control terms have been used. In [21], Nicaise and 
Pignotti examined (jl.ip with g = 0, ao > 0, a\ > and r(t) = r be a constant delay in the case 
of mixed homogeneous Dirichlet-Neumann boundary conditions, under a geometric condition 
on the Neumann part of the boundary. Assuming that < tii < ao, a stabilization result 
is given, by using a suitable observability estimate and inequalities obtained from Carleman 
estimates for the wave equation due to Lasiecka et al. in [14] , However, for the opposite 
case ai > ao, they were able to construct a sequence of delays for which the corresponding 
solution is unstable. The same results were obtained for the case when both the damping 
and the delay act on the boundary, see also [2] for the treatment of this problem in more 
general abstract form. Kirane and Said-Houari [13] considered (jl.ip with ao > 0, a\ > and 
r(t) = t be a constant delay in the case of initial and Dirichlet boundary conditions. They 
established general energy decay results under the condition that < a\ < ao- 

Recently, the stability of PDEs with time- varying delays was studied in [51 [TT 1 [23" 1 [24" 1 [25] . 
In [23], Nicaise et al. analyzed the exponential stability of the heat and wave equations with 
time- varying boundary delay in one space dimension, under the condition < a\ < \/l — dao, 
where d is a constant such that r'(t) < d < 1, V t > 0. In [23J, Nicaise and Pignotti studied 
the stabilization problem by interior damping of the wave equation with internal time- varying 
delay feedback and obtained exponential stability estimates by introducing suitable Lyapunov 



u a -k Au+ / div[a(x)g(t - t)Vu(t)]6.t + b(x)h(u t ) + f(u) = 0, 
Jo 
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functionals, under the condition < y/1 — dao in which the positivity of the coefficient a± 
is not necessary. 

Motivatied by these results, we investigate in this paper system (jl.ip under suitable 
assumptions and prove a general decay result from which the exponential and polynomial 
types of decay are only special cases. Our main novel contribution is an extension of previous 
results from |13[ |2"T] to time- varying delays with not necessarily positive coefficient a\ of the 
delay term. This extension is not straightforward due to the loss of translation-invariance. 
For our purpose, we introduce new energy and Lyapunov functionals, which take into account 
the dependence of the delay with respect to time. 

The paper is organized as follows. In Section [2] we present some assumptions and state 
the main result. The general decay result is proved in Sections El 

2 Preliminaries and main result 

In this section, we present some assumptions and state the main result. We use the standard 
Lebesgue space L 2 (Vt) and the Sobolev space Hq{Q) with their usual scalar products and 
norms. Throughout this paper, Cj is used to denote a generic positive constant. 
For the relaxation function g, we assume the following (see |18t I19j): 

(Gl) g(t) : [0,oo) —> (0,oo) is a non-increasing C 1 function such that 

/•oo 

1 - / g(s)ds = l>0. 
Jo 

(G2) There exists a positive non-increasing differentiable function such that 

g'(t) < -£(t)g(t), t > 0, 

and 

r+00 

/ £(t)dt = 00. 
Jo 

For the time- varying delay, we assume as in [23] that there exist positive constants to, r 
such that 

< r < r(t) < r, V t > 0. (2.1) 
Moreover, we assume that the speed of the delay satisfies 

r'(t) < d < 1, V t > 0, (2.2) 

that 

r G iy 2 '°°([o,r]), VT>0 (2.3) 
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and that ao, «i satisfy 



|ai| < vl — dao- 
As in [23], let us introduce the function 

z(x, p, t) = u t (x, t — r(t)p), x G $7, /)G (0,1), £ > 0. 
Then, problem is equivalent to 



(2.4) 



(2.5) 



uu(x, t) — Au(x, t) + / g(t — s)Au(x, s) ds 
Jo 

+aoUt(x, t) + aiz(x, l,t) = 0, (x, t) G £1 x (0, oo) 

T(t)z t (x, p, t) + (1 - T'(t)p)z p (x, p, t) = 0, 

u(x,t) = 0, fx.f) G <9S1 x fO.oo). ( 2 - 6 ) 

z(x,0,i) = u t (x,t), 

u(x,0) = u {x), u t (x,0) = tii(x), 

, ^(x,p,0) = / (x,-pr(0)), 

We now state, without a proof, a well-posedness result, which can be established by 
combining the arguments of [12\ [13] . 

Lemma 2.1 Let (|2,ip - (|2.4p 6e satisfied and g satisfies (Gl). T/ien given no G Hq(Q), 
u\ G L 2 (0), /o G L 2 (Q x (0, 1)) and T > 0, t/iere exists a unique weak solution («, z) o/ i/ie 
problem (I2.6P on (0, T) suc/i £/iai 

uG C7(0,r;flJ(n))nC 1 (0,T;L 2 (n)), «t G L 2 (0,T;i7 1 (fi))nL 2 ((0,T) x O). 

Inspired by |19} [2"5] , we define the new energy functional as 



(x,p,t) G x (0, 1) x (0, oo), 
(x,i) G dil x [0,oo), 
(x,t) G $7 x (0,oo), 
x G 0, 

(x,p) G O x (0,1), 



E(t) :-- 



u\ + [\ - J g(s)ds) \Vu 



dx + -(goVu)(t) 



4, 



= A(s-i),,2 



u t (x, s)dxds, 



(2.7) 



where £, A are suitable positive constants, and 



(oo W )(t) = 
We will fix £ such that 

2a 

and 



s(t - s)|v(t) - v(s)| 2 dsdx, Vw£ L 2 (ft). 



n jo 



|Ql| 



£ > and £ 



ai 



>o, 



r 



log- 



l a i| 



In fact, the existence of such a constant £ is guaranteed by the assumption (|2,4p . 
Our main result is the following. 



(2. 



(2.9) 
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Theorem 2.2 Let (|2.ip - (|2.4p 6e satisfied and g satisfies (Gl) and (G2). TTten i/iere exisi 
two positive constants K,k such that, for any solution of problem (jl.ip . the energy satisfies 

E(t) < Ke~ kItt o as)ds , V t > t . (2.10) 

Remark 1 Note that the exponential or polynomial decay estimate is only a particular 
case of ()2.10p , More precisely, we can obtain exponential decay for = a and polynomial 
decay for = a(l + t) , where a > is a constant. 

Remark 2 Estimate (|2.10p is also true for t £ [0,io] &2/ virtue of the continuity and 
boundedness of E{t) and£(t). 

3 General decay of the solution 

As mentioned earlier, the proof of the general decay result is given in this section. We have 
the following lemmas. 

Lemma 3.1 Let ()2.ip - (|2,4p be satisfied and g satisfies (Gl). Then for all regular solution 
of problem (jl.ip . the energy functional defined by (|2.7p is non-increasing and satisfies 

E'{t)<hg' oVu){t)-\g{t) f \Vu\ 2 dx - d [ [u 2 t (x, t) + u 2 (x, t - r(t))] dx 
z 1 Jn Jn 

-y/ I e xis - t) u 2 t (x,s)dxds<0 (3.1) 



t-r(t) Jn 



for some positive constant C\ . 

Proof. Differentiating (|2.7|) and by (jl.l|) . we obtain 

E'(t) = u t utt + (l - J g(s)ds\ Vu ■ Vu t - ~g(t)\Vu\ 

+ [ g(t-s) [ Vu t {t) ■ [Vu{t) - Vu(s)]dxds + ]- [ g'{t-s) [ \V u{t) - Vu(s)\ 2 dxds 
Jo Jn 2 ,/ J n 

+ | I u 2 (x,t)dx-| f e- x ^u 2 (x,t-r(t))(l-r(t))dx 

r-t 



dx 

t 



-a| / I e- x(t -^u 2 t (x,s)dxds 
2 Jt-r(t) Jn 

/ ututt + Vn • Vut — / g(t — s)Vu(s) ■ Vut(t)ds 
Jn l Jo 

-\g(t) [ |Vn| 2 dx + i( 5 'oVn)(t) + i I u 2 {x,t)d. 



dx 



| / e- XT ^u 2 (x,t-r(t)){\-r'(t))dx-\^ f [ e- x ^u 2 (x,s)dxds, 
2 Jn 2 J t - T (t) Jn 
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and then, using integration by parts, the assumptions (|2.ip - (|2.2p and some manipulations as 
in [23], 

E'(t) = — oo / u 2 (x,t)dx — at j Ut(t) j ut(t — r(t))dx — —g(t) f |Vu| 2 dx 
Jn Jn Jn 2 

+ l {g i oVu)( t) + i f u 2 (x,t)dx-^ [ e- x ^uUx,t-T(tm-T'(t))dx 
-a| / / e- x ^u 2 t (x,s)dxds 



t-r(t) Jn 

< — ao j u 2 (x,t)dx — a\ I u t (t) I u t (t — r(t))dx — —g(t) j |Vu| 2 dx 
Jn Jn Jn 2 

+ hg' oV«)(t) + | [ u 2 t (x,t)dx-Ul-d)e- x ~ f u 2 (x, t - r(t))dx 



a| /" / e- x ^u 2 t (x,s)dxdi 



2 



t— r(i) Jf2 



<\{g' o Vu)(t) - hg(t) \Vu\ 2 dx - (a 



u 2 (x, t)dx 



2 

-A~ / / e- A(t ~ s) u 2 (x,s)dxds. (3.2) 
2 Jt-r(t) Jn 

Combining (jH%]) - (pl9"]) . (13H and hypothese (Gl), (J3H]) is established. □ 

Now we are going to construct a Lyapunov functional L equivalent to E. For this purpose, 
we define the following functionals: 

I(t) := / uutdx, (3.3) 
Jn 

K(t):=-fu t f g{t - s)(u(t) - u{s))dsdx, (3.4) 
Jn Jo 

Set 

L(t) = NE(t) + eI{t) + K{t) (3.5) 



where N and e are suitable positive constants to be determined later. Similar as in [19] . we 
can prove that, for e small enough while N large enough, there exist two positive constants 
/Si, fa such that 

PiE(t) < L{t) < f3 2 E(t), V t > 0. (3.6) 
The following estimates hold true. 

Lemma 3.2 Under the assumption (Gl), the functional I satisfies, along the solution, 
the estimate 

I\t)<- 1 - [ \Vu\ 2 dx + C 2 f {u 2 {x,t) + u 2 (x,t-T{t))}dx + Cz{goVu)(t). (3.7) 
2 Jn Jn 



G 



Proof. Differentiating and integrating by parts 

l'(t)= / uf dx + / u I Au — / g(t — s)Au(s)ds — aoUt(t) — aiut(t — r(t)) \ dx 
Jn Jn \ Jo J 

= [ u 2 t dx-l I |Vu| 2 dx + f Vu • j g(t - s)(Vu(s) - Vu(t))dsdx 
Jn Jn Jn Jo 

— ao u{t)ut{t)dx — a\ I u(t)ut(t — r(t))dx. (3.8) 
Jn Jn 

Now, using Young's inequality and (Gl), we obtain (see [19] ) 

j Vu- / g(t - s){Vu(s) - Vu{t))dsdx 
Jn Jo 



g(t - s)\Vu(s) - Vu(t)\ds) dx 



2 



f 1 — I 

<5 \Vu\ 2 dx + ——(goVu)(t), V 5 > 0. (3.9) 
Jn 44 

Also, using Young's and Poincare's inequalities gives 

-oq / u(t)u t (t)dx < S I \Vu\ 2 dx + C(5) [ u 2 t dx, (3.10) 
Jn Jn Jn 



-ai I u(t)u t (t - r(t))dx < 5 / \Vu\ 2 dx + C(5) / uf(t - T(t))dx. (3.11) 
Jn Jn Jn 

Combining (|3.8p - (|3.1ip and choosing S small enough, we obtain (|3.7p . □ 

Lemma 3.3 Under the assumption (Gl), the functional K satisfies, along the solution, 
the estimate 

K'(t)<-(! g(s)ds-25] I u 2 t dx + 5 I \Vu\ 2 dx + ^-{g o Vu)(t) 
\Jo J Jn Jn o 

- -4{g'° V«)(t) + 6 I u 2 t {t-T(t))dx. (3.12) 



n 
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Proof. By exploiting (j 1 . 1 j) and integrating by parts, we have 



K'(t) = (l- J g(s)ds\ j ^ Vu- J g{t- s)(Vu(t) -Vu(s))dsdx 

+ f ( [ g(t- s)(Vu(s) -Vu(t))ds) dx- f u t f g'(t - s)(u(t) - u(s))dsdx 
Jn \Jo J Jn Jo 

g(s)ds / u 2 dx + / I / g(t — s)(u(t) — u(s))ds\ \aQUt(t) + a\Ut(t — T(t))]dx. 
Jn Jn \Jo J 



o 

Using Young's and Poincare's inequalities, we obtain (see [191 H] 

"t \ r ft 



I - [ g{s)ds\ [ Vu - [ g(t - s)(Vu(t) - Vu(s))dsdx 
Jo J Jn Jo 

f C 
<8 / \Vu\ 2 dx + — (go Vu)(t), 

Jn o 



f u t [ g'(t-s)(u(t)-u(s))dsdx<5 [ u 2 dx - ^r(g' oVu)(t), 
n Jo Jn ° 



g(t - s)(u(t) - u(s))ds ) [a ut(t) + aiu t (t - r(t))]dx 



n \Jo 



<ir(g oVu)(t) + 5 I u 2 dx + 5 I u 2 t {t -r(t))dx. 



n 
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Combining all above estimates, f)3. 12j) is established. □ 
Now, we are ready to prove the general decay result. 
Proof of Theorem 12. 2L Since g is positive, we have, for any to > 0, 

"to 



ft fto 

/ g(s)ds > / g(s)ds := g > 0, t> t . 
Jo Jo 



By using (|3.ip . (|3.5|) . (\3.7h and (13. 12H . a series of computations yields, for t > to, 

N N f f 

L'(t)<-(g'oVu)(t)--g(t) / |Vu| 2 dx - N& / [u 2 (x, t) + u 2 (x, t - r(t))]dx 
1 1 Jn Jn 



2 

d 
~2 



e x( - s ~ t) u 2 {x,s)dxds + eC 2 / [u 2 t {x, t) + u 2 t (x, t - r(t))]dx 
t-r(t) Jn Jn 

[ \Vu\ 2 dx + eC 3 (go Vu){t) - ( I g{s)ds - 2s) f u 2 t dx + 5 I \Vu\ 2 dx 
Jn \Jo J Jn Jn 



+ Vu)(t) -^f(9'° Vu)(t) + 5 j^u 2 t {t- r{t))dx 

- [(NCi + go) -25- eC 2 ] u 2 dx + (eC 3 + ±C?j (g o Vu)(t) 

+ (f-f)(^v„,«)-(f-*)/ n |v^ 



(NCi - 5 - eC 2 ) I u 2 (x,t-T(t))dx 
Jn 



^ s -^u 2 {x,s)dxds. (3.13) 

't-r(t) Jn 

At this point, we choose e small enough such that e < ^j- and (|3.6p hold, and 5 sufficiently 
small such that 

ai = — - 5 > 0. 

As long as e and 5 are fixed, we choose N large enough such that 

AT (~i 

NCi - 25 > 0, a 2 = Nd - 5 - eC 2 > and a 3 = ^ > 0. 

2 5 

Thus, it follows from (G2) and (l3TT3]l that 

L'(t) <— / u 2 dx — a\ \ \Vu\ 2 dx — a 2 / uf (x, t — r(t))dx 
2 Jn Jn Jn 



2 



t-r(t) Jn 



e K»-t) u 2( Xj s )dxds + a 4 {g o Vu){t) 



< -C 6 E(t) +a 4 {go Vu)(t), V t > t , (3.14) 
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where 04 = eC^ + ^ 



\d > 0. It follows from (j3TT3]h (Gl) and ([H]) that 



£(*)£'(*) < - C 6 S(t)E(t) + 04$ (t)(s o Vu)(t) 

< - c 6 £(t)E(t) - °Vu)(t) 

< - C 6 t(t)E(t) - C 7 E'(t), V t > t , 



(3.15) 



That is 



F ; (t) < -C 8 mE(t) < -k£(t)F(t), V t > t 



(3.16) 



where 



F(t) = £(t)L(t) + C 7 E(t) 



is clearly equivalent to E(t) and A: is a positive constant. 

Consequently, (|2.10p can be obtained by (|3.6p and (|3.16p . □ 
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